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We have obtained the exact ground state wave functions of the Anderson-Hubbard model for 
different electron fillings on a 4x4 lattice with periodic boundary conditions. When compared 
to the uncorrelated ground states (Hubbard interaction set to zero) we have found evidence of 
very effective screening, producing smaller charge inhomogeneities due to the Hubbard interaction, 
particularly at 1/2 filling, and have successfully modelled these local charge densities with non- 
interacting electrons that experience a static screening of the impurity potentials. Further, we have 
compared such wave functions to self-consistent real-space unrestricted Hartree-Fock solutions and 
have found that these approximate ground state wave functions are very successful at reproducing 
the local charge densities, and may indicate the role of dipolar backflow in producing a novel metallic 
state in two dimensions. 



PACS numbers: 71.10.-w,71.23.-k,71.27.-|-a 

Many transition metal oxides display phenomena that 
are believed to be associated with both strong electronic 
correlations and disorder. Various metal-to-nonmetal 
transitions,™ and some properties of the underdoped 
high-Tc cuprate superconductors^i^i^ are examples of 
such physics. Presently the study of such systems is a 
very active field of research in condensed matter physics. 

Theoretically, the simplest model that hopefully rep- 
resents some of the key physics of such systems is the 
so-called Anderson-Hubbard Hamiltonian. The Ander- 
son model^ is given by 

Ha = n,,^ - t (cl^ c^- ^ + h.c)j (1) 

where i,j ~ 1 ... TV denote the sites of the lattice, {i,j) 
implies that i and j are near neighbours, Ci^a i'n-i,^) is 
the destruction (number) operator for an electron at site 
i with spin a, and the hopping energy is denoted by t. 
The on-site energy at site i is given by Vi, and in this 
report we have examined a 50/50 binary alloy model, 
where Vi is set equal to W/2 (for an A site) or —W/2 
(for a B site). The particular complexion of disorder 
that is used in this study can be seen in Fig. 2 - the 
filled circles represent the A sites, while the open circles 
represent the B sites. The electron interactions that are 
included are represented by the Hubbard term, given by 

Hh = u ^^T^^a (2) 

i 

The Anderson-Hubbard model is formed from the sum of 
Ha and Hh- 

There are two dimensionless energy scales in this prob- 
lem: W/t and U/t. Due to the challenge of the computa- 
tional task (see below), in this report we present results 
for U/t = 8 only - such an interaction energy is equal to 
the bandwidth of the non-interacting ordered systems. 



that being 8t. We have examined the weak to strong dis- 
order regimes using W/t = 4,6,7,8,9,10 and 12, and 
have completed calculations for fixed electronic densities 
of 1/4, 3/8 and 1/2 filling. 

In our numerical work we have studied a binary al- 
loy for a 4x4 lattice with periodic boundary conditions, 
and have included an on-site Hubbard repulsion energy 
equal to the bandwidth of the ordered, non-interacting 
problem. Due to the lack of translational periodicity, 
one must determine all of the 16-choose-8 squared (for 
1/2 filling), or roughly 166 million probability ampli- 
tudes, and to date there are no exact wave functions 
available for such a large spatially disordered system. 
(While work on ordered systems^ has indeed studied 
larger Hilbert spaces, previous work'' on this and related 
(spinless fermion) spatially inhomogeneous Hamiltonians 
have not surpassed Hilbert spaces of dimension ~ 3 x 10^, 
less than 1/500*'' the size that we have used.) The role 
of exact diagonalization computational work in the elu- 
cidation of the physics of strongly correlated electronic 
systems is well known - it provides numerically exact 
solutions that are free from approximations. The only 
drawback is that relatively small systems must be exam- 
ined. Larger systems can indeed be studied using quan- 
tum Monte Carlo simulations, and have been completed 
in both two* and three^. dimensions. Our data is com- 
plementary to such work. Therefore, these results can 
now be used as benchmarks gauging the success of pro- 
posed theories, such as Hartree-Fock, which is discussed 
below, or the popularly used dynamical mean field the- 
ory (DMFT).^° Some analogous benchmarking compar- 
isons of DMFT for ordered systems have recently been 
completed. 

For this cluster one can exactly diagonalizc the Hamil- 
tonian matrix and find the ground state wave function 
using the Lanczos algorithm. Besides the energy of these 
states, we have found the local charge densities as well 
as one characterization of the magnetic properties given 
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TABLE L In (a) the exact local charge densities are listed 
for the U/t — non- interacting ground state at 1/2 filling for 
W/t = 4; (b) shows analogous data, but now for the exact 
interacting systems with U/t = 8; (c) shows the local charge 
densities for the interacting problem, now calculated within 
the Hartree-Fock approximation. 

by Si • Sj for i and j being nearest neighbours. 

The disorder potential without the Hubbard energy is 
sufficient to produce a strongly inhomogeneous system, 
even for moderately small disorder. This is demonstrated 
in part (a) of TABLE |T] in which we have listed the local 
charge densities for W/t = 4 and J7 = 0. The A sites 
have charge densities from 0.14 to 0.28, whereas the B 
sites have local charge densities from 1.71 to 1.84. How 
effective is the Hubbard interaction energy in screening 
the disorder potential? For U/t — 8 the answer from 
exact diagonalization is found in part (b) of TABLE [J 
With the Hubbard interaction the local charge densities 
for both A and B sites now range from only 0.91 to 1.11, 
and inspection of the data for the non-interacting and 
interacting charge densities makes evident the effective- 
ness of the Hubbard energy in producing a much more 

I 



where the effective local fields hf are given by 

ht^iS^) K^{S~) ■ (5) 

Then, one must find self consistently the local spin- 
resolved charge densities and local fields that minimize 
the variational estimate of the ground state energy. 



uniform charge density. 

However, it is to be emphasized that when one is away 
from 1/2 filling the homogenization of the local charge 
densities is far from complete. That is, making the Hub- 
bard energy of the order of or larger than the disorder 
potential does not necessarily lead to complete screening 
of the disorder. As an example note that for C//t = 8 and 
W/t = 4 for 3/8 filling, one finds a bimodal distribution 
of the local charge densities peaked at roughly 0.5 (for 
the A sites) and 1.0 (for the B sites). 

One important element of work towards the under- 
standing of such systems is the characterization of elec- 
tronic screening of the (impurity) disorder potential, and 
recent publicationsi^ii^ have stressed the importance of 
disorder screening in interacting systems. We have ex- 
amined the simplest variant of such an approach. To be 
specific, for non-interacting electrons we have determined 
the ground states with screened impurity potentials, viz. 
for effective on-site energies given by 

V, - V^ff = ^ (3) 

Solving for the local charge densities for such a model, 
we fit the static dielectric constant, e, for each filling 
and W/t by minimizing the mean-squared difference of 
these densities in comparison to the exact local charge 
densities. This procedure is quite successful, especially 
(as expected) at 1/2 filhng, and in TABLE |TT] we com- 
pare charge densities for W/t=12 from which one can see 
the impressive agreement (note that this is for U/t = 8) 
between these approaches. For this parameter set and 
filhng, we find e « 2.6. 

The simplest and most common approach to finding 
approximate solutions for interacting systems is that of 
using a self-consistent, real-space unrestricted Hartree- 
Fock (HF) approach, which previously have been em- 
ployed in various ways in previous studies^ i^^'^^ . 

In such a HF approach one ignores terms that are pro- 
portional to fiuctuations about mean values squared, and 
thus one approximates the local Hubbard interaction be- 
ing replaced by (here we ignore the possibility of local 
superconducting pairing correlations) 



(4) 

I 

The HF approach is a variational one in which the 
lowest energy solution within the product state Hilbert 
space is found. A comparison of the energies for different 
fillings and disorder is shown in Fig. [U and the agreement 
between the exact and HF energies is seen to improve 
with increasing disorder, as expected. 
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TABLE 11: In (a) the exact local charge densities are given for 
the U/t = 8 interacting ground state at 1/2 filling for W/t = 
12; (b) shows analogous data, but now for the exact non- 
interacting system in which the impurity potential W/t —12 
is screened as in Eq. ([3}, for a static, homogeneous dielectric 
constant of £ = 2.6. 



It is common rubric that HF works best for disordered 
systems (in comparison to ordered systems), but the cri- 
tiquing of such a belief requires that one be in possession 
of exact solutions with which one can compare approx- 
imate solutions. One such comparison is found in part 
(c) of TABLE U shows the local charge densities for the 
HF ground state for W/t = 4 and U/t = 8 for 1/2 fill- 
ing (the spins are found to be coUinear for this state, 
so it is permissible to set hi = - see below). Clearly 
the agreement between the exact and HF ground states 
is very good - the maximum absolute difference between 
exact and HF charge densities is 0.019. This quantitative 
agreement is found to be repeated for all W/t, U/t and 
electronic fillings. One particular example of the success 
of this approximation is found for W/t = 12, C//t = 8 at 
1 /4 filling - for these parameters the largest absolute dif- 
ference between the exact and HF local charge densities 
is only 0.002. 

We now turn to what one can learn from the HF so- 
lutions concerning the spin degrees of freedom for such 
systems, both the spin-spin correlations and the spatial 
arrangement and orientations of the spins. For the HF 
solutions there is some freedom of choice in specifying 
the possible effective fields that are allowed. That is, if 
the hf are set to zero then the only nonzero expectation 
values of (S^ • Ca), where a — x,y,z, are for the z com- 
ponents. Such solutions correspond to (potentially) non- 
paramagnetic spin arrangements with coUinear spins. If 
the hf are allowed to be non-zero but constrained to be 
real then — and both the x and z components 
can be nonzero. This allows for the possibility of the HF 
solutions having non-coUinear (but still coplanar) spins 
lying in the xz plane - below we will refer to such so- 
lutions as "twisted spin configurations". Both of these 
circumstances lead to restricted HF solutions; only when 
the hf are allowed to assume complex values does one 
have the fully unrestricted solutions. In our numerical 
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FIG. 1: [Colour online] The exact vs. real-space self- 
consistent Hartree-Fock energies for the binary alloy Hamilto- 
nian discussed in the text. The Hubbard interaction is fixed 
to he U/t = 8, and the lower/middle/upper curves correspond 
to fixed electronic densities of one-quarter /three-eighths/half 
filling, namely 8, 12 and 16 electrons (8e, 12e and 16e). 



work on a 4x4 lattice, we have found that for all sys- 
tem parameters and fillings the HF states correspond to 
either coUinear or coplanar spin configurations. 

In contrast to the success of the HF wave functions in 
producing an excellent approximation for the local charge 
densities of the exact ground state, since these wave func- 
tions are product states the spin degrees of freedom are 
essentially classical (although the moments on each site 
can be different than 1/2) and therefore the effectiveness 
of this approximation in reproducing the spin correlations 
in the ground state is not expected to be as good. Indeed 
this is what we find - there is only a reasonable qualita- 
tive similarity between the exact and HF near-neighbour 
spin-spin correlations. As one typical example note that 
for W/t = 6 and U/t — 8 at 1/2 filling, the average (ab- 
solute) difference between the exact and HF evaluations 
of {Si ■ Sj) for i,j being near neighbours is almost 20% 
- for some near neighbours the HF correlations do not 
even have the same sign as found in the exact solutions. 
Further, comparing the exact and HF results using a spa- 
tially averaged mean-squared difference, while the local 
charge densities improve with increasing disorder, sim- 
ilar to the energies of Fig. [TJ the spin-spin correlations 
become increasingly worse with increasing disorder. 

The spin degrees of freedom that are present in the HF 
solutions do contain some potentially important informa- 
tion, and here we mention one possible relationship of our 
results to other studies. That is, the spin arrangement 
that we find may have a bearing on a recent proposal for 
a novel metallic phase in two dimensions The authors 
of this paper found that at 1/2 filling and for a disorder 
strength that was approximately equal to the Hubbard 
energy there was a significant delocalization of the ef- 
fective one-electron HF eigenstates at the Fermi energy. 
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FIG. 2; The spin configuration for a 1/2 filled system with 
W/t = 7 and U/t = 8. The filled circles denote A sites and 
the open circles denote B sites. The length of each vector is 
proportional to the magnitude of the magnetic moment on 
that site. 



and thereby forecast the existence of a metallic state in 
two dimensions. An understanding of the physics of this 
result remains an outstanding problem. 

We have found an interesting result in our HF studies 
that relates to this problem, that being that when we 
focus on 1 /2 filling, only for W/t =7 and 8 (for U/t = 8) 
do we find non-coUinear spin configurations; for all other 
values of W/t we find coUinear HF ground states. In 
Fig.[2]we show the non-coUincar spin configuration found 
for 1/2 filling for W/t = 7 and U/t = 8. It is in this 
same parameter regime that these authors^'^ found the 
existence of a novel metallic phase, and the role that these 
spin configurations may play in such a transition is an 
interesting question to explore. Further, some evidence 



of this physics is also found in a comparison of the exact 
and HF wave functions. Out of all of the exact wave 
functions that we obtained, the only (non-singlet) triplet 
ground state was for 1/2-filling and W/t — 7. In fact, for 
this system the HF state had a net magnetic moment of 
1, consistent with the exact ground state. 

To answer the question of what, if anything, does the 
non-coUinear spin configuration of the HF solution have 
to do with the enhanced delocalization of electrons found 
in Ref. 15, we draw a parallel to the work studying a sin- 
gle hole moving in an antiferromagnet background (de- 
scribed by the t — J model). Speculation of the potential 
physics behind the metallization in two dimensions found 
in Ref. 15 can then be gleaned from the dipolar back- 
flow analysis of a vacancy moving in an antiferromagnetic 
background. That is, the spin "twists" are generated by 
allowing for the mobile vacancy to minimize its energy. 
Such a state produces a non-zero quasiparticle weight, 
itself a requirement for a metallic phase. 

Summarizing our results, we have completed the exact 
diagonalization for a 4x4 cluster with periodic bound- 
ary conditions for a binary alloy with a repulsive on-site 
Hubbard interaction. Our analysis has shown that the 
interactions generate an effective screening of the disor- 
der potential for some fillings and disorder. Further, our 
real-space self-consistent Hartree-Fock (HF) solutions are 
seen to be able to reproduce quantitatively the local 
charge densities with considerable success, but not the 
quantitative details of the spin correlations. Our HF so- 
lutions display non-coUinear, coplanar spin twists in the 
ground state at 1/2 filling when the disorder and interac- 
tion energies are comparable. This report will be followed 
by a longer paper with many more numerical results and 
detailed analysis. 
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